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■^ ' Abstract 

ri ' We investigate solid-solid and solid-liquid transitions of elastic flexible off- 

lattice polymers with Lennard- Jones monomer-monomer interaction and 
'"^3 ■ anharmonic springs by means of sophisticated variants of multicanoni- 

£^ ' cal Monte Carlo methods. We find that the low-temperature behavior 

O . depends strongly and non-monotonically on the system size and exhibits 

broad similarities to unbound atomic clusters. Particular emphasis is ded- 
icated to the classification of icosahedral and non-icosahedral low-energy 
polymer morphologies. 
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O 1 Introduction 

O' 

The understanding of the crystallization behavior of single flexible and elastic 
polymer chains is of particular interest in natural and technological molecu- 
lar building-block systems, where functionalization is based on the individual 
structural properties of small molecules forming larger-scale composites. Among 
d ' these are, e.g., hybrid systems of solid substrates and "soft" molecules, where 

it is known that the formation of regular structures can be mediated by the in- 
teraction with the amorphous or crystalline surface of the substrate [1] . Tuning 
the self-assembly of individual polymers near substrates is relevant for molecu- 
lar nanoelectronic applications and materials design which has become a huge 
field of applied research. As an example from life sciences, capsomer proteins 
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of spherical viruses assemble to form a perfectly icosahedral capsid hull, which 
shelters the viral genome [2]. Thus, frozen structures of single polymer chains 
can serve as the basic elements of larger assemblies on nanoscopic scales. As we 
will show here, the crystallization behavior of elastic polymers exhibits strong 
similarities to the cluster formation of colloidal (or atomic) particles. 

The general behavior of a polymer in solvent has already been subject of 
numerous studies [3, 4]. At high temperatures or in good solvent the interactions 
of different parts of the polymer arc of little relevance and due to entropic 
effects the dominating structures are extended random coils. This changes at 
the tricritical 0-point where the polymer behaves in three dimensions effektively 
like a Gaussian chain, up to logarithmic corrections [5]. The effective repulsion 
caused by excluded-volume effects and the monomer- monomer attraction are 
in perfect balance. In the globular regime, the chain is collapsed and compact 
conformations dominate which, having almost no internal structure, resemble 
a liquid state. At lower temperature, energy optimization leads to a freezing 
towards crystalline structures. In the thermodynamic limit, this transition is 
expected to be of first order, whereas for the O-collapsc a second-order phase 
transition is predicted [6, 7]. 

Under certain conditions, the two transitions also may coincide in the ther- 
modynamic limit as it was shown in a recent investigation of a bond-fluctuation 
model [8]. However, this is not a generic feature of flexible [9] and elastic poly- 
mers [10]. 

While the coil-globule transition is well analyzed and further progress in the 
analysis of logarithmic corrections to the scaling behavior would require the 
investigation of extremely long chains, less is known for the freezing transition 
especially in the case of off-lattice polymers. Due to the intrinsic similarities of 
the applied model to colloidal clusters it is useful to interpret the crystallization 
of elastic polymers with reference to atomic Lcnnard- Jones clusters and to use 
the corresponding nomenclature. In fact, if one replaces 'atoms' by 'monomers' 
many statements from the cluster field remain also true for polymers. Besides 
the numerous studies that have been performed on the freezing transition and 
ground-state properties of Lennard- Jones clusters, there are only a few publi- 
cations on Lennard- Jones homopolymers. The authors of Ref. [11] were able to 
show principal relationships between both classes of systems and motivated the 
transfer of terms and concepts from the cluster case to the field of polymers. 
Recent studies [12, 13, 14, 15, 16] mainly focused on the general behavior of 
elastic polymers but did not refer in detail to the polymer length dependence of 
the freezing transition, the ground-state conformations or solid-solid transitions. 

The aim of this work is to close this gap and to compare the low-temperature 
properties of polymers and clusters. The rest of the paper is organized as 
follows. In the subsequent section we will introduce the model and mention 
some problems that arise in the simulations. We then give a short review of the 
thermodynamic properties of Lennard- Jones clusters in Sec. Ill, followed by the 
introduction of a geometrical order parameter in Sec. IV. After the description 
of the method and the explanation of a few technical details in Sec. V, we will 
present and discuss our results in Sec. VI. Finally, the main findings of the paper 



are briefly summarized in the concluding part. 



2 Model 

We employ a model for elastic flexible polymers, where all monomers interact 
pairwisely via a truncated and shifted Lennard- Jones (LJ) potential: 

t£3° d (r«) = FLj(min( Wc )) - V L] (r c ). (1) 

Here, r^ denotes the distance between the ith and jth monomer, r c is the cutoff 
distance, and 

F LJ (r)=4e[(a/r) 12 -(a/r) 6 ] (2) 

the standard LJ potential. In the following, e is set to I and a = ro/2 1 / 6 with 
the minimum-potential distance r$ = 0.7. In this work we used and compared 
three different cutoff distances r c = 2.5a, 5c, oo. 

Covalent bonds are modeled by an additional FENE-potential [17] for adja- 
cent monomers: 

V FEnE (rii+i) = -yi? 2 ln{l - [(r«+i - r )/R} 2 }. (3) 

The potential possesses a minimum coinciding with tq and diverges for r — > 
ro ± R with R = 0.3 here, if is a spring constant set to 40. 
Eventually, the polymer energy is given by: 

N itf-1 

E = 2 £ VlTW + £ ^FENE(r M+1 ). (4) 

3 A Short Review of Lennard-Jones Clusters 

Much work has been done in LJ cluster studies to identify and classify global 
energy minima. At least for N < 200 the ground states are known with 
high certitude [18, 19, 20]. Most global minima are constructed via a com- 
mon scheme: usually an icosahedral core is covered by an incomplete overlaycr, 
which gets more populated with increasing system size N until a magic number 
(N = 13,55,147,309,561,923 [21]) is reached and the next complete icosahe- 
dron is formed. Two types of overlayers occur: If only a few monomers are 
added to the core they gather at one of the faces filling the available space most 
efficiently. This leads to an hep- like or anti-Mackay packing [Fig. 1(a)]. With 
increasing system size, other faces are occupied and additional atoms adsorb at 
the corners. Thereby atoms on different faces interact only weakly. Hence, a 
different type of structure is energetically favored when a certain system size is 
exceeded. The external monomers adapt the partly fcc-like structure of the core 
and build a Mackay overlaycr, which also allows the occupation of the core edges 
at the expense of non-optimal distances [Fig. 1(b)]. With a few exceptions, one 




Figure 1: (Color) Overlayer (bright) types on an icosahedron's face (dark): (a) 
anti-Mackay (hep), (b) Mackay (fee) [18]. 



finds anti-Mackay ground states for clusters with iV =14-30,56-81,85, whereas 
sizes N =31-54,82-84,86-146 lead to Mackay global energy minimum (GEM) 
[18]. For some systems, whose sizes differ notably from the optimal values 
N = 13,55,147,..., non-icosahedral geometries become energetically competi- 
tive. The ground state of the 38-atom cluster is an fee-like truncated octahedron 
[22], systems of 75-77 [23] and 102-104 [24] atoms form Marks decahedra with 
a central fivefold symmetry axis and the cluster with 98 atom forms a Leary 
tetrahedron [25]. 

When considering non-zero temperatures one finds that below the melting 
temperature clusters occupying Mackay-like conformations undergo a solid-solid 
transition where the overlayer transforms to anti-Mackay type [26] . Thereby the 
transition temperature increases with system size. Clusters with non-icosahedral 
GEM typically change at very low temperatures to icosahedral structures with 
larger population at intermediate energies. 

4 Geometrical Classification of Solid Elastic Poly- 
mers 

Structural changes between different geometries appear to be a characteristic 
feature of the cluster's behavior at low temperatures. Since we expect a similar 
behavior also for clastic polymers, we are interested in classifying the geomet- 
rical states of a conformation in order to specify conformational transitions. 
While many studies refer to bond orientation parameters [27] that are calcu- 
lated via spherical harmonics, we restrict ourselves to informations provided by 
the contact map which can be updated instantly during the simulation without 
further computational expenses. We consider two monomers as being in contact 
if their distance is smaller than a threshold r contac t. Due to this definition also 
bonded monomers do not need to be in contact. The total number of monomer 
contacts is not an appropriate measure for classification since in the interior of 
a frozen polymer every monomer has usually exactly 12 neighbors. Instead, the 
contacts between that 12 neighbors reflect their arrangement which corresponds 
to the local geometry. In Fig. 2(a-d), different conformations of a monomer and 
its 12 neighbors are shown. Counting the contacts between the neighbors is a 
simple but efficient way to characterize different types: Only the icosahedral 
cell [Fig. 2(a)] reveals 30 shell contacts corresponding to the 30 edges of an 
icosahedron. It always appears in the center of an icosahedral conformation. 
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Figure 2: (Color) Possible shell conformations (r c = 2.5a) of a monomer pos- 
sessing 12 (a-d) or 11 (e) neighbors, (a) icosahedron, (b) elongated pentagonal 
pyramid [28], (c) cuboctahedron (fee), (d) triangular orthobicupola [28] (hep), 
(e) incomplete icosahedron. Sticks illustrate shell contacts, not bonds. 



Consequently, if there is no such basic element, the global geometry cannot be 
icosahcdral! On the other hand, icosahedral cells are also formed by a suffi- 
ciently large anti-Mackay overlayer at the corners of the icosahedral core. If the 
number of outer monomers is too small, one might find instead the defected 
icosahedral cell [Fig. 2(e)], a monomer with 11 neighbors forming 25 shell con- 
tacts. The total number of both structures nj c is a suitable "order" parameter 
which allows a classification of the global geometry at low temperatures, given 
roughly by 

= non-icosahedral 

= 1 icosahedral + Mackay, (5) 

> 2 icosahedral + anti-Mackay. 

More precisely, if ni c = 0, the polymer forms a non-icosahedral structure, e.g., 
it is decahedral or fcc-like; tt.j c = 1 indicates icosahedral geometry with Mackay 
overlayer or a complete icosahedron which might possess a few monomers bound 
in anti-Mackay type. Finally, for n- lc > 2, the monomers form an icosahedral core 
with a considerably extended anti-Mackay overlayer. The probabilities p nic (T) 
for the different values of n[ c as a function of temperature provide the necessary 
information to reveal structural transitions. 

Figure 2(b) shows the elongated pentagonal pyramid with 25 shell contacts 
which is the basic module of five-fold symmetry axes in icosahedra and deca- 
hedra. It also occurs along the edges of an icosahedral core, which is covered 
by an anti-Mackay overlayer. Hence, it appears in icosahedral conformations 
with TV > 31 and decahedral structures. Besides, it is formed at the edges of 
the central tetrahedron in conformations with a tetrahedral symmetry. An ex- 
ample is the ground state of the cluster with N = 98 or low-energy minima for 
N = 159 and N = 234. In consequence, the total number of elongated pen- 
tagonal pyramids n cpp can be used to distinguish decahedral and tetrahedral 
structures. Figures 2(c) and (d) show conformations which hardly differ because 
both of them possess 24 neighbor-neighbor contacts. They occur in almost all 
geometries considered here. Only the truncated fec-octahedron (i.e., the ground 
state of the 38mer) does not exhibit triangular orthobicupolae [Fig. 2(d)]. Nev- 
ertheless, since cuboctahedra [Fig. 2(c)] are related to fee- and triangular or- 
thobicupolae to hep-packing, their observation is of some use in other cases. In 
this study, we will mainly focus on the analysis of the number of complete and 
defected icosahedral cells. n w . 
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Figure 3: Density of states for systems forming complete icosahedra. 



5 Simulation Method 



For the precise simulation of this model we have to face several challenges. 
Firstly, since the elastic bonds represent an additional restriction which causes 
multiple energy barriers, it is necessary to generate an appropriate simulation 
dynamics that allows the tunneling through these barriers and fast conforma- 
tional changes in order to reduce autocorrelation times. We achieve the latter 
by applying three different conformational updates. The first causes a rela- 
tively small change by a simple shift of a monomer, the second propagates a 
monomer along the chain to a new position which is defined within a suitable 
local coordinate system, whereas the third alters only the linkage by swapping 
two bonds between four nearby monomers. Details of these update procedures 
will be reported elsewhere [29]. 

Secondly, investigating the thermodynamic behavior of the polymer within 
the temperature interval of interest requires the sampling of regions of energy 
space with spectral densities differing by many hundreds to thousands orders of 
magnitude (Fig. 3). In order to enable the system to perform a random walk 
in energy space we applied the multicanonical method [30], where in iterative 
runs the multicanonical weight function W{E) is estimated [31]. In a final 
production run a conformation X contributes with a generalized Boltzmann 
weight e~ /3E ^W(E('X.)) such that all energies are roughly equally probable. 
The density of states can then be determined by 
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where H(E) is the histogram of the production run. Thereby, the Metropolis 
[32] acceptance probability for a change from system conformation X to X' is 



replaced by 

w* x '> = °»° ('■ e-" <Em - B>xn} wim$) ■ (7) 

The parameter /3 influences the simulation only at the beginning of the weight 
estimation whereas it is canceled out when the weight function converges. To 
determine the multicanonical weights we performed several hundred single sim- 
ulations, each with 10 5 N updates. The final production run consisted of at 
least 10 9 N updates. The obtained histogram can be reweightcd to any desired 
temperature as long as the corresponding energy interval is covered by the sim- 
ulation data, yielding the canonical energy distribution for this temperature. 

Finally, the shape of the energy landscape might prevent the exploration of 
a narrow global energy minimum with standard Monte Carlo methods in cer- 
tain cases - a problem which is well-known from LJ clusters consisting of 38, 
75-77, 98 or 102-104 atoms [33, 34] - so that further improvements of the sim- 
ulation method are required. As described in the previous section, the number 
of icosahedral cells rii c (X) allows for the geometrical classification of conforma- 
tions X. Based on n; c (X), parameters ^(X) can be introduced to divide the 
conformational space in different parts. If non-icosahedral minima shall be ex- 
plored, conformations without icosahedral cells must be treated differently from 
icosahedral conformations. Using 

KX) " 1 2 if n ic (X) > 1, (8) 

the statistical weight of the two fractions can be enhanced or suppressed within 
the framework of multicanonical Monte Carlo sampling by employing multiple 
weight functions. Within an extended multicanonical ensemble a conformation 
X is than represented with a probability 

p(X) oc e-^ x >T^ (X )(£(X)) (9) 

that depends on -E'(X) and v(X.) such that (7) changes to 

-faccept(X, X ) = 

-/3(B(X')-B(X)) W v( x')(E(X')) 



W„w(E(X)))- 

Consequently, different geometries which correspond to different values of v can 
be tuned to participate equally at any energy and free-energy valleys of different 
depths and widths can be sampled easily. Therefore, the investigation of the 
solid-solid transitions is no longer a problem of barriers and huge autocorrelation 
times. While following the branches of different geometry, the only remaining 
task is to reach energies that correspond to a temperature which is sufficiently 
smaller than the temperature of the solid-solid transition. 



In a similar way the system can be enabled to change more frequently be- 
tween Mackay and anti-Mackay conformations. By using the parameter 

/,-v-s / 1 if «ic(X) < 1 , , 

"&> = { 2 if n ic (X) > 2 ( U ) 

and introducing corresponding weight functions, the simulations can be tuned 
to hit both types of structures with equal feqeuency at any energy as long as 
such conformations exist. 



6 Discussion and Results 

6.1 Ground states 

Global energy minimum conformations of the investigated polymers reveal wide 
similarities to ground-state configurations of LJ clusters, i.e., for almost all sys- 
tem sizes, the ground state is of icosahedral type (Fig. 4). At zero temperature, 
the monomer- monomer bonds cause only small deviations since the minimum 
of the FENE potential is close to the equilibrium distance. Besides, the chain 
can arrange within the mostly unaltered monomer conformation to minimize 
the bond potential. As a result, bonds between different shells of an icosahe- 
dral core are rare since the corresponding monomer distances are smaller than 
the equilibrium distance and entail higher bond energies. For longer chains, 
the central icosahedral cell is strongly compressed and one end of the chain is 
usually located in the center of the cell, thus avoiding the inclusion of the most 
inappropriate distance (from the central monomer to one of its neighbors) in 
the chain twice. A similar effect occurs in decahedral conformations. Bonds 
between monomers on the central axis are favorable since their length is close 
to the optimum. This forces the polymer chain to adapt this axis at low tem- 
peratures [Fig. 5(a), (c)]. 

We find perfect icosahedra for system sizes N = 13,55,147,309. If the 
system size exceeds this 'magic' sizes, the polymer builds an icosahedral core 
with an anti-Mackay ovcrlayer which grows with the chain length. At some 
point (N > 30 and N > 80), the over layer adopts the structure of the core by 
changing to the Mackay type. Further increase of N completes the outer shell 
and leads to the next icosahedron. A few polymers of certain sizes show ground- 
state conformations that correspond to different non-icosahedral geometries. We 
find a truncated octahedron for N = 38 and a decahedral configuration for 
N = 75 — 77 [Fig. 5(a)]. Some deviations are caused by the cutoff of the LJ 
potential, so the chains with N = 81,85,87,98, 102 monomers possess lowest- 
energy structures that do not correspond to the respective clusters unless an 
untruncated LJ potential is applied. Using the cutoff r c = 2.5a we find for 
N = 81,85,98,102 icosahedral ground states with Mackay ovcrlayer and for 
N = 87 a conformation with two merged icosahedral cores (Fig. 6). 



N=13 



N=55 




N=147 






Figure 4: (Color) Icosahedral (putative) ground-state conformations for different 
system sizes. 




Figure 5: (Color) Non-icosahedral ground-state conformations: (a) TV = 75, 
r c = 2.5cr decahedral, (b) N = 98, r c = 5er tetrahedral, (c) N = 102, r c = 5<r 
decahedral. 




Figure 6: (Color) The structure of minimal energy for TV = 87, r c = 2.5a is 
formed by two entangled icosahedral cores. The ends of the chain coincide with 
the respective centres. 



6.2 Thermodynamics of complete icosahedra 

The four investigated chains that form complete icosahedra (N = 13, 55, 147, 309) 
exhibit a very clear and uniform thermodynamic behavior. We observe two sep- 
arate conformational transitions indicated by peaks in the specific C(T) heat 
and the fluctuations of the radius of gyration d(r gyr )/dT. As shown in Figs. 7 
and 8, the icosahedra melt in the interval 0.3 < T < 0.5 and a liquid-like regime 
is reached where the polymer arranges still in a globular shape but exhibits no 
distinct structure. Hence, the icosahcdral order parameter {n- lc ) s=s 1 changes 
to (?ii c ) ~ 0. The corresponding peak in the normalized specific heat increases 
rapidly with system size and allows in principle a precise determination of the 
melting temperature. However, all solid-solid or liquid-solid transitions consid- 
ered here must not be understood as thermodynamic transitions in a strict sense 
since all investigated systems are small and dominated by finite-size effects. For 
longer chains one would expect, in analogy to the thermodynamic behavior of LJ 
clusters, the crossover from icosahedral ground states to decahedral (N > 1500) 
and later to fcc-like [35] structures (N > 20000) which exhibit a different crys- 
tallization behavior. Therefore the extrapolation to the thermodynamic limit 
by means of finite-size scaling is unfortunately not an appropriate choice. It is 
worth noting that for N = 309 the result differs considerably from the pure LJ 
cluster in which case a recent study [36] found two clearly separated peaks. 

Increasing T further leads to the collapse transition at temperatures 1 < 
T < 2. Above this temperature the chains arrange randomly in extended con- 
formations. For almost all investigated systems this cross-over does not produce 
a clear signal in the specific heat except for the chain with N = 309 where a 
shoulder emerges at T ~ 1.8. If the chain is shorter this energetic signal of the 
collapse transition is suppressed by finite-size effects. However, geometric quan- 
tities like the radius of gyration and its fluctuation (Fig. 8) give some insight. 
It is obvious that the solid-liquid transition remains well separated from the 
coil-globule collapse, moreover the intermediate temperature interval increases 
within the investigated chain length interval. 

The solid phase is dominated by the extremely stable icosahedron which, 
however, depending on the temperature, exhibits surface defects. Although the 
mobility of a great majority of the monomers is strongly restricted, there are 
still changes which can be observed. Firstly, the chain, i.e., the linkage of the 
monomers still transforms relatively freely. Only at extreme low temperatures 
the lowest energy (chain-) conformation gains thermodynamic relevance. Start- 
ing from the center, the number of bonds connecting different shells is reduced. 
At T = the first three shells are connected only by a single bond, so that for 
N = 13, 55, 147 one end of the chain is the central monomer and the other is 
at the surface. For N = 309 bonds between the corners of the third and the 
fourth shell exist also at T = 0, since the length of this bonds is roughly equal 
to alternative bonds within the third shell, which include the corner monomers. 
Secondly the entire icosahedron undergoes a compactification with decreasing 
temperature. 
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Figure 7: The specific heats for chains forming complete icosahedra (r c = 2.5a). 
Peaks indicate the melting transition at 0.3 < T < 0.5, whereas the coil-globule 
collapse shown in the inset is less pronounced. 

6.3 Liquid-solid transitions of elastic polymers 

In the following, we analyze the behavior of small elastic polymers in the liquid- 
solid transition regime. Particular emphasis will be dedicated to the chain- 
length dependence of geometric changes in the conformations of the polymers 
while passing the transition line. 

The size dependence of the specific heat for polymers with sizes up to N = 55 
is illustrated in Fig. 9(a). The ground state of the short 13mer is the energet- 
ically stable icosahedron whose almost perfect symmetry is only slightly dis- 
turbed by the FENE bond potential. The melting transition is indicated by a 
the peak at T = 0.33. With increasing system size the low-temperature con- 
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Figure 8: The normalized fluctuations of radii of gyration for chains forming 
complete icosahedra give clear evidence for the freezing and the collapse transi- 
tion as well as their increasing separation. 
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Figure 9: Specific heat for (a) 13 < TV < 55 and (b) 55 < N < 147, r c 
For better visibility, the curves are shifted by a constant offset. 



2.5a. 



formations are less symmetric and the peak becomes broader. For N > 30, 
the crossover from Mackay-like ground states to anti-Mackay conformations 
takes place, the corresponding peak increases with growing N and shifts from 
T ~ 0.042 to higher temperatures. Finally, for N = 55, there is only one melt- 
ing transition left and anti-Mackay-like structures are strongly suppressed. As 
an effect of the bond elasticity, the melting transitions occur at slightly higher 
temperatures than in the case of pure L J clusters [26] . 

The behavior of chains containing between 55 and 147 monomers [Fig. 9(b)] 
again generally corresponds to that of LJ clusters [26] with the differences that 
the melting transition occurs at higher temperatures and that the Mackay-anti- 
Mackay transition temperature increases much faster with system size. Besides 
there are a few chains with ground states of types different to L J clusters which is 
induced by the truncation of the LJ potential as mentioned above. We find anti- 
Mackay ground states up to system sizes of N = 80. In contrast to the clusters, 
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Figure 10: Whereas the Mackay-anti-Mackay crossover is not recognizable in 
the specific heat curves, the probability of occurrence of specified numbers of 
icosahedral cells, p„ ic , reveals the transition temperature (r c = 2.5a). 



N = 81 and N — 85 possess a GEM of Mackay type. Hence we encounter a 
solid-solid transition also in these cases. We do not find peaks of the solid-solid 
transition in the specific heat for TV = 81,82,85, and the crossover to anti- 
Mackay conformations can only be identified in changes of the structural "order" 
parameter n ic defined in (5). In Fig. 10 we therefore also show the probability 
p„ ic that rii c < 1 or n- lc > 2 as a function of temperature. While one would 
expect higher transition temperatures for growing system size, this prediction 
fails in the case of N = 85, where the anti-Mackay energy minimum is almost 
as deep as the GEM due to an optimal arrangement of the outer monomers. 
Nevertheless, for N > 90, the transition shifts rapidly to higher temperatures 
(Fig. 11) and manifests in the specific heat as well. Whereas in the case of 
pure LJ clusters, the two peaks remain separated up to sizes of 130 atoms 
[26], we observe both transitions merging already for N ps 100. In contrast 
to the polymers, the Mackay-anti-Mackay transition temperature of clusters 
even decreases near N = 120, presumably because the anti-Mackay overlayer 
is almost complete, leading to a spherical and therefore stable conformation. 
This indicates that in this polymer model anti-Mackay structures loose weight 
compared to atomic clusters. For N = 75, one finds the crossover between 
decahedral and icosahedral conformations indicated by a small peak at T w 0.8. 
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Figure 11: Already for N > 110 the entire solid phase is dominated by Mackay 
(n; c = 1) conformations (r c = 2.5a). 
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Figure 12: (Color) Conformations of minimal energy for N = 85 with (a) anti- 
Mackay and (b) Mackay over layer. 

6.4 The cutoff's influence 

At this point we cannot yet judge whether the differences to the behavior of 
atomic clusters at medium temperatures are caused by the truncation of the LJ 
potential or by the polymer topology, i.e., the bonds. To answer this question, 
we will discuss three interesting cases with modified cutoff. First, we consider the 
polymer TV = 85 with the original LJ potential, i.e., r c = oo. The corresponding 
cluster is the largest in the interval 55 < N < 147 with an anti-Mackay ground 
state [18]. Since we use the unaltered LJ potential, any deviation to the cluster 
behavior is caused by the bond potential only. In contrast to the previous 
simulations with truncated potential, we retain the anti-Mackay ground state 
where the 30 outer monomers completely cover 10 faces of the icosahedral core 
and build an energetically favored structure [Fig. 12(a)]. One might notice that 
there are no bonds connecting monomers on different faces directly since the 
bond length would be too far from the potential minimum. This means that 
for very low energies only a few bond configurations are allowed and that the 
anti-Mackay state is much less metastable than the Mackay state [Fig. 12(b)] for 
which many more bond configurations with low energies are possible. This leads 
to an entropic dominance of the latter in the temperature interval 0.002 < T < 
0.08 as visible in Fig. 13, a behavior that has never been reported for atomic 
clusters. In spite of large statistical errors we find a signal in the specific heat 
at T « 0.002, since at this temperature the energetically favored anti-Mackay 
state prevails. The transition back to anti-Mackay conformations extends over 
a wider temperature interval around T « 0.08 and cannot be localized in the 
specific heat. Again the change from a LJ cluster to a polymer leads to a greater 
prominence of Mackay conformations. 

To prove the influence of the cutoff on the thermodynamics at higher temper- 
atures we doubled the cutoff to r c = 5a rj 3.12, minimizing differences to the un- 
truncated potential and investigated the two systems N = 98, 102. For the cor- 
responding clusters, two separate transitions where observed at medium temper- 
atures which stands in contrast to the polymer with original cutoff (r c = 2.5a) 
where both transitions merge. Applying the enlarged cutoff we observe only a 
slight shift in the peak positions whereas the major differences to the LJ clusters' 
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Figure 13: The polymer with N = 85 and untruncatcd LJ potential (r c = oo) 
changes from anti-Mackay (T < 0.002) to Mackay conformations (0.002 < T < 
0.1) and back (0.1 < T < 0.34). 

behavior as demonstrated in Ref. [26] persist (Figs. 14, 15). It is still impossible 
to distinguish both conformational transitions clearly by means of the specific 
heat, i.e., the temperature domain where anti-Mackay like conformations prevail 
is very small. We may conclude that this difference to cluster behavior (i.e., the 
suppression of anti-Mackay states) is mainly an effect of the bonds. 

The second interesting outcome of the simulations are the different ground 
states. As in the case of unbounded clusters with untruncated interaction, we 
obtain a GEM of tetrahedral geometry for N = 98 [Fig. 5(b)] and a decahe- 
dral ground state for N = 102 [Fig. 5(c)]. We can also identify the solid-solid 
transition to Mackay structures: For N = 98, (n; c ) ps changes as expected 
to (tMc) « 1 at T « 0.0075 (see inset of Fig. 14) and a tiny bulge is formed 
which, however, has to be considered with caution, since statistical errors in- 
crease strongly at such low temperatures. In the other case (TV = 102) the 
decahedral-icosahedral crossover occurs at a higher temperature T w 0.02 (see 
Fig. 15), with a more prominent signal in the specific heat. 

After all we see some evidence that during the change from atomic LJ clus- 
ters to polymers, Mackay-like structures are more favored independently of a 
truncation of the L J potential. This results in a shift of the Mackay-anti-Mackay 
(or Mackay overlayer melting-) transition to higher temperatures and the for- 
mation of a Mackay dominated temperature interval for the 85mer if the original 
LJ potential is employed. 

A further clear effect of the truncation of the potential arises in the case 
of TV = 87. This particular number of monomers allows the formation of two 
merged icosahcdral cores, each with a icosahedral center of high density (Fig. 6). 
The overall shape of this conformation reminds of a cylinder and is therefore not 
optimal for the original long-range LJ potential which tends to form spherical 
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Figure 14: (Color) Although there is almost no change in the specific heat the 
values p„ ic indicate clearly the crossover to a non-icosahedral ground state if the 
cutoff r c is increased (N = 98). 



s 





Figure 15: (Color) Same as Fig. 14 for N = 102. The solid-solid transition to 
the decahedral ground state leaves a small but distinct peak in the specific heat. 
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Figure 16: For N = 87 the transition to a unique GEM of excentric shape is 
visible in the specific heat, in p„ ic as well as in the radius of gyration r gyr which 
strongly increases during the transition. 




Figure 17: (Color) Low-energy conformations of the polymer with N = 234 
(r c = 5er) with (a) icosahedral, (b) tetrahedral, and (c) decahedral geometry. 



structures. Applying the cutoff r c = 2.5er, the short-range interactions gain 
importance and the double-core conformation appears to be the ground state. 
At T as 0.12 the system undergoes the cross-over to anti-Mackay structures 
and also Mackay conformations play a minor role (Fig. 16). Since thereby the 
polymer becomes more spherical, we observe a decrease of the radius of gyration. 
Note that in the ground-state conformation (Fig. 6) both ends of the chain are 
located in the centers of the two cores, thereby minimizing the chain energy in 
the way mentioned above. 

6.5 Non-icosahedral minima of larger systems 

To demonstrate the capability of our method to find different structures of 
minimal energy we now present results for a relatively long chain. Using a 
5cr-cutoff, we investigated the 234mer which has the suitable length to form 
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tetrahedral as well as decahedral minima. To distinguish the two cases, we used 
the number of elongated pentagonal pyramids [Fig. 2(b)] n epp in addition to 
rii c for the classification. While n cpp < 6 for conformations with decahedral 
symmetry, n epp = 24 in the tetrahedral case (ri; c = for both). We were 
able to find the expected structures for each geometry (Fig. 17) and a final 
minimization with the conjugate gradient technique [37] reveals that while the 
icosahcdral minimum represents the GEM the energy difference to the minima 
of the two other geometries is very small (see Table 1). This is also the case 
when considering the minimized respective cluster conformations, i.e., neglecting 
bonds and applying the untruncated LJ potential. The obtained value for the 
icosahcdral LJ cluster minimum is in perfect agreement with an earlier study 
[19], hence we may assume that also the energies of the decahedral and the 
tetrahedral state are close to the true minimal values. 



Table 1 : Minimal energy values for conformations of different geometry of the 
polymer with N = 234 (r c = 5a) and the 234-atom LJ cluster. 

icosahcdral tetrahedral decahedral 

polymer -1458.877 -1458.431 -1458.036 

LJ cluster -1465.924 -1465.529 -1465.116 



7 summary 

In this study, we demonstrated the capability of the multicanonical sampling 
method to investigate the complete thermodynamic behavior of flexible elastic 
polymer chains containing more than 300 monomers within reasonable time on a 
standard workstation. In order to achieve this we developed new conformational 
updates which enhanced the performance significantly. Analyzing the behavior 
of the radii of gyration it became evident that the solid-liquid transition and 
the coil-globule collapse will remain well separated also for much longer chains. 
In our study of conformational properties we focused on the low-temperature 
regime. It turned out that the extensive results of the last decades' research 
on the subject of atomic Lennard- Jones clusters provide an excellent framework 
for also understanding the liquid-solid transition behavior of the polymers. To 
observe the expected low-temperature transitions including conformations of 
completely different geometry, we introduced the total number of icosahcdral 
cells nic as an order parameter for the classification and distinction of the re- 
spective structures. Fortunately, this new quantity did not only prove useful in 
visualizing the conformational changes, but also enabled us to explore the tetra- 
hedral and decahedral minima for TV = 75, 98, 102, and 234. For this purpose we 
introduced multiple weight functions which is a relative simple modification to 
the standard multicanonical method which should be in principle also adaptable 
to other methods like parallel tempering. 

The low-temperature behavior of elastic polymers exhibits strong similar- 
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ities to atomic Lcnnard- Jones clusters. In almost all cases, the global en- 
ergy minimum possesses an icosahedral core and an overlayer of anti-Mackay 
(hep) or Mackay (fee) type. Differences in the putative ground states (e.g., for 
N = 81,87,98, 102) could be ascribed to the truncation of the Lennard- Jones 
potential. Furthermore, in our analysis of the long-range cutoff influence on 
crystallization we made manifest for N = 87 and r c = 2.5a the solid-solid tran- 
sition from icosahedral double-core to single-core conformations and for N = 75, 
r c = 2.5cr and N = 98, 102, r c = 5a the cross-over from non-icosahedral to icosa- 
hedral structures. As a common effect of the bonds, all melting transitions were 
shifted to higher temperatures compared to atomic Lennard- Jones clusters. Be- 
sides, for chains of length 55 < N < 147 further deviations arose: While in 
the case of clusters for 82 < N < 130 the Mackay-anti-Mackay transition was 
clearly separated from the core-melting this only holds for polymer chains with 
81 < N < 100. For longer chains both transitions merge. For the chain with 
N = 85 and untruncated Lennard- Jones potential, we find a temperature range 
where Mackay conformations dominate, embedded in regions of anti-Mackay 
structures. We regard this last two findings as evidence for the tendency of the 
bond potential to suppress anti-Mackay structures. 
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